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ABSTRACT 

We show that comparisons of He II Lya forest lines of sight to nearby quasar populations can strongly 
constrain the lifetimes and emission geometry of quasars. By comparing the He II and H I Lya 
forests along a particular line of sight, one can trace fluctuations in the hardness of the radiation field 
(which are driven by fluctuations in the He II ionization rate). Because this high-energy background 
is highly variable - thanks to the rarity of the bright quasars that dominate it and the relatively short 
attenuation lengths of these photons - it is straightforward to associate features in the radiation field 
with their source quasars. Here we quantify how finite lifetimes and beamed emission geometries affect 
these expectations. Finite lifetimes induce a time delay that displaces the observed radiation peak 
relative to the quasar. For beamed emission, geometry dictates that sources invisible to the observer 
can still create a peak in the radiation field. We show that both these models produce substantial 
populations of "bare" peaks (without an associated quasar) for reasonable parameter values (lifetimes 
~ 10 6 -10 8 yr and beaming angles < 90°). A comparison to existing quasar surveys along two He II 
Lya forest lines of sight rules out isotropic emission and infinite lifetime at high confidence; they can 
be accommodated either by moderate beaming or lifetimes ~ 10 7 -10 8 yr. We also show that the 
distribution of radial displacements between peaks and their quasars can unambiguously distinguish 
these two models, although larger statistical samples are needed. 

Subject headings: cosmology: theory - intergalactic medium - quasars: absorption lines - quasars: 
general 



1. INTRODUCTION 

Quasar spectra provide powerful probes of both the 
source's properties and the intergalactic medium (IGM). 
One particularly useful aspect is the so-called "proxim- 
ity effect," which describes the highly-ionized zone sur- 
rounding each bright quasar. Measuring the transition 
from this zone to the more uniform ionizing background 
characteristic of the average IGM provides an estimate of 
the magnitude of that background as well as source prop- 
erties like the quasar luminosity, lifetime, variability, and 
emission geometry. 

The classical proximity effect test uses the zone along 
the line of sight to a bright quasar by measuring the 
H I Lya forest. The increased radiation background 
causes excess transmission in the forest near to the 
quasar, which on e can model to extract the ionizatio n 
rate of H I, r H i (|Baitlik et all 119881: IScott et al.ll2000D . 
However, the H I proximity zone typically only spans 
a few Mpc, within the overdense neighborhood of the 
quasar's massive host galaxy, and it is difficult to dis- 
entangle the intrinsic transm ission bias of this region 
(jFaucher-Giguere et al.ll2008D . 

Another flavor, the so-called "transverse proximity ef- 
fect," can be more powerful. Here, we measure the im- 
pact of radiation from a foreground quasar on a dif- 
ferent Lya forest skewer: as the line of sight passes 
the foreground quasar, the forest should show an en- 
hanced ionizing background. However, if the foreground 
quasar's light is beamed - toward the observer, in this 
case - or if the quasar has a short enough lifetime, 
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its radiation may not even intersect the Lya forest 
skewer at all, or at least it may not strike it at the 
point of closest approach. Thus the pattern of en- 
hancements near foreground quasars should reveal in- 
formation abo ut the lifetime, beaming, and variabil- 
ity of quasars (ICrottsI Il989t iMoller k Kiaergaardlll992t 
lAdelbergerll2004t ISchirber et alJl2004t ICrofdl200¥ . 3 

With the H I Lya forest, one still only expects to see 
an enhancement in the biased e nvironment very close 
to the quasar (within < 3 Mpc) (jFaucher-Giguere et al.1 
120081 iHennawi et all 12006). Such close pairs of quasars 
are rare, even in modern data sets. Only recently have 
large surveys and more careful analysis led to the de- 
tection of the transverse effect: iGoncalves et al.l (2008) 
used higher-ionization metal line ratios to measure life- 
times ~ 3 x 10 7 yr with no evidence for anisotropic emis- 
sion, while K irkman k Tvtlen (|2008| ) searched 130 quasar 
pairs, with separations < 3 Mpc, and found no increased 
transmission in the foreground but decreased transmis- 
sion in the background, which they explain by appealing 
to rapid (~ 10 6 yr) variability. 

In this paper, we argue that the He II Lya forest 
offers a better measurement of the transverse proxim- 
ity effect. It has two important advantages. First, 
the He II ionization rate (for which we will use T as 
a shorthand for rH e u) is much more variable than Thi 
dFardal et aTJIl 998: Ma selli k Ferrara 2005; B olton et ail 
2006; Mciksin 20Q3 iFurianettB [200930). Each quasar 
dominates the local high-energy radiation field in a much 
larger region - tens of comoving Mpc - than it does for 
H I, so the signature is much easier to identify. Second, 

3 The classical proximity effect is less sensitive to these effects 
because all of the information comes from a single skewer (and 
hence a single light path). 
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we can compare the He II and H I forests in order to 
measure the hardness of the ionizing background. This 
eliminates the dependence on density (and temperature), 
which affect both the He II and H I fractions in the same 
way, and r emoves any ambiguity due to the quasar en- 
vironment (jFaucher-Giguere et al .1120081 ). To the extent 
that the H I-ionizing background is uniform, we can also 
therefore cleanly measure the He II-ionizing background. 

These two factors make the detection and extraction 
of physical quantities much simpler with the He II forest. 
Indeed, the transverse proximity effect h as already been 
detec t ed along two different lines of sight (j jakobsc n et al.l 
12001 IWorseck fc Wisotzkil 12001 IWorseck et al.l 120071) . 
even with low signal-to-noise spectra. Here we will 
show that the method provides great discriminating 
power between quasar models with finite lifetimes and/or 
anisotropic emission. 

This paper is organized as follows. In Sections [2] and 
[3] we use toy models of quasars with finite lifetimes and 
beaming to compute the fraction of peaks in the radiation 
field that have nearby observable quasars. In Section @] 
we show that the angular distribution of the peak-quasar 
associations can distinguish these two sets of models. Fi- 
nally, we compare to existing observations in Section [5] 
and conclude in Section [6] 

All distances are quoted in proper units unless other- 
wise specified. 

2. QUASARS AND RADIATION PEAKS: FINITE LIFETIMES 

Here we will use toy models and simple statistics to 
illustrate the power of the He II transverse proximity ef- 
fect. Our basic approach is to measure the fraction of 
peaks in the radiation field identified along a Lya forest 
line of sight with observable quasars sufficiently nearby 
to be clearly identified as the peak's source. For isotrop- 
ically emitting quasars with infinite lifetimes, all peaks 
must have such an association, but some will disappear 
if quasars are beamed (because many of the sources will 
not be oriented toward the observer) or if they have finite 
lifetimes (because many of the sources will have shut off 
already). 

For the purposes of this paper, we take a very sim- 
ple toy model in which only a single (nearby) quasar 
induces fluctuations in the He II-ionizing background; 
we assume that the cumulative background from other 
sources is uniform . This is certainly a simplification 
(|Furlanettol [2009a), but we will defer a more compre- 
hensive model to future, m ore detailed work. Figure 1 of 
iFurlanetto fc Dixonl (|2010D shows that this is a reason- 
able approximation in most cases of interest (i.e., near 
strong peaks in the radiation field). 

We define a "peak" in the radiation field as any point at 
which the He II-ionizing background intensity J > a (J), 
where a is a constant. We will usually choose cr = 1, 
because the median ionizing ba ckground is a factor ~ 3 
less than the mean (|Furlanettoll2009aI) . 

In practice, such peaks will be identified as minima 
in the "hardness ratio" 77 = Na e u/N}n cx Thi/T, the ra- 
tios of column densities in singly-ionized helium and neu- 
tral h ydrogen within a given absorber ([Miralda-Escudel 
1993;. We implicitly assume that Thi is uniform, so 
77 cx 1/r oc 1/J. T his appears to be an exc e llent approx - 
imation at z < 4 (jMeiksin fc White! [200l ICroftl [200l . 
This approach has the added benefit of eliminating any 



uncertainty in the density or temperature structure of 
the forest. Note, however, that regions so close to their 
sources so as to be within both the proximity zones of 
He II and H I may be missed by this technique. 

We then define r p as the maximum transverse distance 
within which a quasar can sit and still produce a radia- 
tion peak. If the quasars are Poisson-distributed and the 
mean free path of an ionizing p hoton is rp, the averag e 
radiation field is (J) = 3JV J* (|Meiksin fc White! 12004 , 
where Nq is the average number of quasars inside one at- 
tenuation zone (with radius ro), J* = (L) /(47rro) 2 , and 
(L) is the mean quasar luminosity (ave raged over the 
lumin osity function, which we take from iHopkins et al.l 
[200l . 

There is unfortunately considerable disagreement on 
the mean free path of these photons. At z ~ 2.5, 
models predi ct tha t rp ~ 45-200 comoving Mp c (e.g., 
iBolton et all 120061 : IFaucher-Giguere et al.l 120091) . One 
method uses the IGM density distribution in simula- 
tions and calibrates to the observed average optical depth 
(|Bolton et al.l l2006). and it repro duces the abundance of 
H I L yman-limit systems from IStorrie-Lombardi et al.l 
(1994). However, it relies o n ad hoc geometric as- 
sumptions about the a b sorbe rs. iFardal et al.1 (|1998l) and 
IFaucher-Giguere et alj (|2009| ) use photoionization mod- 
eling of absorbers in the H I forest; the latter repro- 
duces the more recent mean free path measuremen ts (for 
photons that ioniz e H I) of iProchaska et alj ((2009) and 
iSongaila fc Cowiel (|2010j) . However, this method relies 
sensitively on the H I forest absorbers with TVhi ~ 10 15 - 
10 17 cm -2 , which are very difficult to measure, and othe r 
estimates give much smaller values (Fardal e t al.lll998l) . 
The more recent data is calibrated to z ~ 3.7, so it also 
cannot account for possible evolution to the redshifts of 
interest, z ~ 2.5. Neither method takes into account the 
substantial fluctua tions of the He II ionizing background 
(|Furlanettoll2009afl . 

Other arguments suggest that r lies in the middle of 
this range. First, let us assume that the abundance of 



Lya forest absorbers follows a power law cx N„t , close 
to observational limits (e.g.. IFardal et all 11998 ). In the 
optically thin limit with uniform radiation backgrounds, 
the hardness ratio r\ is spatially constant. We can there- 
fore use the distribution of H I absorbers to relate the 
mean free paths of H I and He II ionizing photons (rni 
and rp respectively): 



rm 



(1) 



IFaucher-Giguere et al.l (|2009l ) estimat e that rm = 
85([l + z]/4)~ 4 Mpc, which matches the IProchaska et all 
(2009) measurement at z = 3.6 very well, although the 
redshift evolution is uncertain. In that case, 77 ~ 40-80 
dShull et al.ll2004IZheng et alj|2004t iFechner et aHl2006t 
iFechner fc Reimersl l2007h implies ro ~ 110-160 Mpc at 
z~2.5. 

A separate argument comes from matching the (mea- 
sured) emissivity of the quasar population with the ob- 
served optical depth of the He II forest. In the fluctuat- 
ing Gunn- Peterson approximation, ro > 100 Mpc over- 
produces the ionization ra te compared to observations 
ijDixon fc Furlanettoir2009l ). although there is an uncer- 
tain correction factor in this model. 

We therefore take ro = 45 comoving Mpc as a fiducial 
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tg (yr) Opening Angle (Degrees) 

Fig. 1. — Fraction of radiation field peaks with direct quasar associations for models with finite lifetimes (a) and beamed emission (b). 
In each panel, the thick solid curve shows our fiducial model (see text), the thick dotted curve takes a higher survey luminosity threshold, 
and the thick dashed curve shows the predictions for a survey spanning only 12 comoving Mpc around the Lyce forest skewer and in which 
only peaks with a > 8 are identified. In (a), the (upper) thin curves take ro = 200 comoving Mpc. (This parameter does not affect the 
results with beaming.) In (b), the thin curves show the absolute fractions of all quasars that are visible to the observer (short-dashed) and 
that are not visible but produce measurable peaks in the forest (dotted). 



model but show results for larger values as well. 4 Then 
we find TV, 



ro 



L/L, where 



3N <7 (L) 



5 x 10 12 (T 



45 Mpc 



(2) 



is a characteristic luminosity. (Here the factor 1 + z 
converts our physical coordinates into comoving units.) 
Quasars must sit inside this region and intersect the 
Lya forest skewer in order to produce an observable 
peak. For reference, a bright quasar (~ 10 12 L Q ) has 
r p ~ '"o/vS ~ 20 comoving Mpc with our fiducial value 
for ro; r p ~ 12 comoving Mpc for ro = 150 comoving 
Mpc. 

If quasars have a finite lifetime tQ, the induced radia- 
tion peak along a nearby line of sight can be displaced 
radially from the quasar location, and - because of light 
travel time delay - the quasar also may not be visible 
when the peak is observed. To compute the probabil- 
ity to have a "bare" peak (i.e., without an associated 
quasar), we note that the quasar and peak are both visi- 
ble 1) once the peak has moved within a distance r < r p 
of the source and 2) before tQ has elapsed. The first 
corresponds to a minimum time after the quasar appears 
of 



tn 



(3) 



where r± is the quasar's impact parameter from the Lya 
forest skewer. Before i m i n , the peak along the Lya forest 
skewer is at least r p in front of the quasar and is invisible. 

4 An alternative approach is to scale the proximity zone to the 
He II ionization rate, T. In that case, we remove the uncertainty 
in ro but replace it with equivalent uncertainty in T. 



Of course, if t min > tQ, the peak never reaches the prox- 
imity zone during the quasar's lifetime, so there cannot 
be an association. 

After the quasar shuts off, the peak will continue to 
drift farther from the observer as the corresponding light 
travel time delay increases. The peak will reach the back 
of the proximity zone after a time delay, relative to the 
last quasar photon seen by the observer, of 



r 2 p -r\ 



(4) 



If all quasars have a fixed luminosity, so that r p is 
constant across the population, then the fraction of peaks 
with quasar associations is Fit = N pq / (N p b+N pq ), where 



N pq = n I dr±r±c[t Q - t min (r±)] and (5) 
Jo 
rr P 

N pb = TT / dr±r±c[At mBX (r±) - max(i min - t Q , 0)Q6) 
Jo 

are respectively the volumes within which quasars can 
sit and have a visible peak (or not). Here r max enforces 
the requirement that i m j n < tQ] it is 



(7) 



if fq = ctq/r p < 1 and unity otherwise. The final fac- 
tor in Npi, accounts for the time lag between the quasar 
shutting off and the peak reaching the proximity zone, 
in cases where tQ < t m - m . 

The integrals in equations J!])-© can be performed 
analytically; the results scale roughly as r p from the 
cylindrical geometry of the problem, with a suppression 
at large luminosity because of the finite lifetime limits 
(even though bright quasars have large proximity zones, 
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they can still shut off before the peak becomes visible). 
Of course, we must actually integrate N p b and N pq over 
the quasar luminosity function. Because these factors 
scale like oc L, the most luminous quasars are by far 
the most important for generating peaks in the radiation 
field. 

Figure [TJa shows the resulting fractions for several dif- 
ferent mock surveys. In order to better mimic real sur- 
veys, we take a minimum luminosity threshold when cal- 
culating the fraction of observed associations, though not 
when calculating the total number of peaks (i.e. we limit 
L in the numerator of F\ t but not in the denominator). 
We do not (yet) limit the spatial extent of the survey, in- 
stead assuming that all quasars can be identified out to 
the appropriate r p (L), which is in principle measurable 
from each quasar's luminosity. In practice, this would 
require a survey that is deepest near the Lya forest line 
of sight. 

The thick solid curve takes L > 10 11 L Q and ro = 45 
comoving Mpc at z — 2.5. The fraction of associations 
increases rapidly at tq < 3 x 10 7 yr and then flattens 
out at larger lifetimes. For small lifetimes, we expect the 
result to be ~ Tq , which is the fraction of the proximity 
zone for which the light travel time is less than tQ . This 
is indeed roughly correct; once tq > 1, the curves flatten 
significantly because the more distant quasars only pro- 
vide observable peaks for brief windows of time anyway. 

The dotted curve shows how the fraction varies with 
the survey depth, taking L > 10 12 L Q . Clearly a shal- 
lower survey strongly reduces the number of observed 
associations. However, note that decreasing the limit be- 
low 10 11 Lq has very little effect. One need only identify 
those quasars responsible for strong peaks, which are pri- 
marily bright and moderate luminosity sources. To pro- 
duce a peak, faint sources must already be so close to the 
line of sight that their available volume is small. How- 
ever, although fewer peaks have associations, the varia- 
tion of the curves with tQ is relatively constant with L, 
so a wide, shallow survey may be just as effective as a 
deeper one, if one is confident enough about modeling 
the fainter quasar population. 

The thin solid curve takes ro = 200 comoving Mpc; we 
find that the fraction of associations is roughly propor- 
tional to Tp cx r^ 1 at short lifetimes. A larger attenua- 
tion length increases the fraction of peaks with quasars, 
because the additional sources illuminating each point ef- 
fectively decrease r p (and hence time delay effects) in or- 
der to overcome the background from the other sources. 
The dependence is substantial, so a more accurate esti- 
mate of the mean free path will be essential for detailed 
constraints. 

Finally, the long-dashed curv es assume a survey com- 
parable to lWorseck et al.l (|2007| ): again the thick and thin 
curves take ro = 45 and 200 comoving Mpc, respectively. 
We take L > 10 11 L Q , only include peaks with a > 8, 
and only identify quasars within 12 comoving Mpc of 
the Lya forest skewer. These factors decrease the de- 
pendence on tQ (and hence make the survey less sensi- 
tive) , largely because of the high peak threshold: quasars 
must be very close to the line of sight in order to produce 
such strong peaks. These nearby sources do not provide 
much constraining power, because the light travel time is 
then small compared to tQ ; surveys for bright quasars at 



distances near r p are most efficient. Nevertheless, short 
lifetimes (tg ~ 10 6 yr) would imply very few observable 
associations with this kind of survey (see Section [5]), and 
very long lifetimes tQ > 10 8 yr) would imply almost per- 
fect association. 

3. QUASARS AND RADIATION PEAKS: BEAMING 

We now switch focus to models with infinite lifetime 
but anisotropic quasar emission (or "beaming" ) . For con- 
creteness, we will use a simple biconical emission model, 
in which two oppositely-directed beams each have open- 
ing angle f2. 

In this case the calculation is conceptually simple: how 
often does a beam that remains invisible to the observer 
intersect the Lya forest line of sight between the front 
and back edges of the proximity zone? Some visible 
quasars will have such intersections, some will have none, 
and some invisible quasars will still cause peaks; our 
accounting must include all these possibilities. (Note 
that quasars whose beams do not intersect either the 
observer's line of sight or the Lya forest skewer remain 
entirely invisible and can be ignored.) 

We use a Monte Carlo model to compute these proba- 
bilities. 5 The thin curves in Figure [lb show (1) the frac- 
tion of quasars that are visible to the observer and pro- 
duce visible peaks (short-dashed curve) and (2) the frac- 
tion that are invisible to the observer but still produce a 
peak (dotted curve). The former simply increases with 
£1, of course. 6 The latter initially increases (as quasars 
become more likely to intersect the skewer) and then de- 
creases (as quasars become more likely to be seen by the 
observer) . 

To generate predictions relevant to observations, we 
must include quasars of all luminosities, as before. How- 
ever, if we assume (as in our fiducial model) that we 
detect all quasars within their respective r p (L), then 
the results become independent of luminosity because 
each quasar is treated identically (this differs from the 
finite lifetime case, where tq introduces a second phys- 
ical scale). If, however, r max is fixed (as in a real sur- 
vey), then the integration over the luminosity function 
becomes necessary. Similarly, if the survey is not in- 
finitely deep, the minimum luminosity threshold intro- 
duces luminosity dependence because some associations 
with faint quasars will be missed. 

The thick solid curve in Figure \T]p shows results for a 
survey with r max = r p (L) and L > 10 11 L Q . The num- 
ber of bare peaks can be substantial; roughly two out of 
three peaks will have no visible association if il ~ 7r/2, 
with the ratio declining very rapidly at smaller opening 
angles. As in the finite lifetime shallower survey 

misses associations. However, in this case the shape of 
the curve as a function of il does change with the lu- 
minosity threshold, becoming slightly more sensitive to 
variations in il as the survey deepens. 

The long-dashed curve again shows a survey similar to 

5 Without the limits imposed by the proximity zone, an analytic 
calculation is straightforward - in fact every visible quasar pro- 
duces a peak somewhere - but requiring some part of the beam to 
strike the skewer within r p of the quasar's location makes such a 
model unwieldy. 

6 Note again that all quasars produce "peaks" somewhere along 
the line of sight, but we do not count them unless they are within 
the proximity zone. 
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Fig. 2. — Cumulative distribution function of the peak angle 9 for quasars with finite lifetimes (left) and biconical beaming (right). In 
the beaming case, the distributions are independent of distance from the quasar. In the finite lifetime case, we show several different impact 
parameters between the Lyo forest skewer and the quasar, scaled to the light travel time over the quasar's lifetime, ctQ. 



IWorseck et al.l (|2007j ). which only includes strong peaks 
and nearby quasars. Because there is no physical scale 
other than r p in the problem, simply increasing a - which 
affects all quasars equally - does not affect the fraction of 
bare peaks. On the other hand, limiting r max does, be- 
cause it eliminates the possibility of finding more distant 
luminous counterparts. However, with the large a im- 
posed here nearly all peaks are sourced by quasars within 
the surveyed region, so the finite area makes only a small 
difference to the final curve, which is mostly determined 
by the depth of the survey. 

It is worth emphasizing that the beaming case is some- 
what more robust to observational uncertainties, because 
it does not depend on ro (at least in our simplified 
model). 

4. PEAK-QUASAR ASSOCIATIONS 

So far we have shown that the relative number counts 
of "bare" peaks and those with associated quasars de- 
pends on the lifetime and beaming angle. Here we briefly 
show that the relative radial locations of the peaks and 
their associated quasars can provide additional evidence 
and in particular distinguish these two scenarios. 

Clearly, if each quasar emits isotropically over an infi- 
nite lifetime, it illuminates every point in the universe, 
so the brightest point along a nearby Lya forest skewer 
(at a radial distance from the observer r pea k) will lie at 
the same radial distance as the quasar: i.e., if 9 is the 
angle between the Lya forest skewer and the ray join- 
ing r P eak and the quasar, then 9 = tt/2 for every source, 
barring complex radiative transfer effects and errors in 
localizing the quasar and peak. Deviations from this sim- 
ple expectation therefore indicate more complex quasar 
properties, such as finite lifetimes or beaming. 

4.1. A Finite Lifetime 

With a finite lifetime, a quasar visible through both its 
direct emission and through its influence on a point in the 



IGM must satisfy a time delay criterion: the difference 
in light travel time along these two paths At must be no 
greater than tQ. In terms of the angle 9, this time delay 



At = (r - d LOS )/c = (1 - cos 9)r/c = 



1 — cos t 
sin 9 



( 8 ) 

where r is the total distance from the quasar to the near- 
est point that it illuminates and cZlos is the radial dis- 
tance between that same point and the quasar. This is an 
increasing function of 9, so points closest to the observer 
have the least delay, 7 and it increases monotonically as 
9 approaches tt/2. 

We will construct the probability distribution of the 
angle 9 in a finite lifetime model (recall that it is a 
delta function at 9 = tt/2 for the fiducial infinite lifetime 
model). If the quasar turned on a time t in the past, 
only points with At < t are illuminated, and the peak 
is at the point with the largest 9. Assuming a uniform 
distribution of quasar ages, the cumulative distribution 
of peak locations is therefore 



P«9) = ^=( 1 -^ 
V J to \ sine? 



CtQ 



The left panel of Figure [2] shows this distribution 
for several impact parameters r±, scaled to ctQ w 
30(£q/10 8 yr) Mpc. Two points are immediately obvi- 
ous. First, 9 < tt/2 always; once the quasar illuminates 
this point of absolute closest approach, the peak must 
remain there until the quasar shuts off. Thus we expect 
a clear asymmetry between the forward and backward 
directions. Second, the distribution depends strongly on 
the impact parameter from the skewer. Nearby quasars 
have short delays, so 6 = tt/2 is most common. But 
those at relatively large impact parameters only rarely 

7 We assume here that the quasar is effectively at an infinite 
distance from the observer. 
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have the timing just right to attain 9 = n/2. Thus faint 
quasars - which must be nearby to influence the radia- 
tion field - provide little additional information on finite 
lifetimes, and surveys for bright quasars over wide areas 
are most productive. 

Note that here we have not restricted ourselves to 
quasars with visible peaks within a distance r p from the 
source. The importance of this restriction depends on 
the impact parameter r±_ ; those peaks with cos 9 ~ 1 
will be difficult to identify in practice. 

4.2. Quasar Beaming 

Now we consider the distribution of 9 in a beaming 
model with infinite lifetime; specifically, biconical emis- 
sion. More complicated emission geometries are of course 
possible and can dramatically change the estimates in 
this section. 

For a quantitative picture, we again turn to the dis- 
tribution of 9. In the beamed case, quasar-peak pairs 
visible to the observer will typically not have 9 = n/2, 
because the beams can only subtend (at most) an angle 
fi from 9 = or n (along the radial direction) . The right 
panel of Figure [2] shows the resulting cumulative prob- 
ability distributions of 9 for a range of opening angles 
(generated with a Monte Carlo model). 8 

As expected, when each beam has Q = n (i.e., isotropic 
emission), the peak is always located perpendicular to 
the line of sight. As ft decreases, the probability of this 
configuration decreases rapidly, and by il = 7r/2, it is 
vanishingly rare. Instead, the peak drifts farther and far- 
ther from the perpendicular, because the illuminated re- 
gion is oriented more and more directly toward (or away 
from) the observer. 

There are two clear differences from the finite life- 
time case. First, the beaming results are distance- 
independent, because there is no physical scale in the 
problem (although "peaks" beyond r p will no longer be 
visible in practice). Second, the distributions are sym- 
metric about 9 = it/2, because we assume two beams 
directed in opposite directions; 9 > n/2 indicates that 
different beams illuminate the skewer and observer. The 
second difference therefore depends on the particular 
model of quasars, but the first will unambiguously dis- 
tinguish finite lifetimes from geometric effects. 

5. COMPARISON TO EXISTING OBSERVATIONS 

So far, searches for neighboring quasars have been 
conducted along two of the five lines of sight with 
He II Lya forest data. The best sample so far sur- 
rounds the line of sight to HE 2347 -4342, whose He II 
forest has been extensively studied (IZheng et al.l 120041: 
Shull et al.ll200l lFechner_fc Refiners! 120071 : IShull et al l 
20101 ) . IWorseck et al.1 (|2007ft searched for nearby quasars 
with L m i n > 1-2 x 10 Lq and a maximum im- 
pact parameter ~ 12 comoving Mpc. They discovered 
two nearby quasars, both with associated peaks in the 
high-energy radiation background (actually identified as 
troughs in the hardness ration ij). 9 Three other regions 

8 Note that we normalize these distributions to unity, so we 
include only those sources that are both visible and create peaks 
along the line of sight; the other configurations discussed in Section 
[3] are not relevant for this test. 

9 They also report a third quasar just below their survey limit, 
which we do not include to ensure completeness. However, because 



of the spectrum have a radiation field at least as hard as 
these (with a — 8), so the fraction of associa tions is 2/5. 

The second skewer is toward Q0302-003. iHeap et all 
(2000) identified a strong tra nsmission feature near z = 
3.05. iJakobsen et al.l (|2003[ ) subsequently identified a 
quasar at a projec ted distance of just 1.77 com oving Mpc 
from this feature. IWorseck fc Wisotzkl (|2006l) later per- 
formed a deeper search (to comparable magnitude and 
volume limits as above) and identified one more quasar 
within ~ 5 Mpc (projected) of the He II Lya forest line 
of sight and a third inside the classical proximity zone 
of Q0302-003 (we ignore the last, since its effects are 
d ifficult to disentangle) . 

IWorseck fc Wisotzkil (|2006l ) find peaks in the hard- 
ness of the radiation background near the locations of 
all of these quasars. They do not quantify the num- 
ber of peaks without quasars, but visual inspection of 
their Fig. 7 show that the sample also contains two other 
hardness peaks of comparable amplitude to those with 
known quasars. If these are included (so that the frac- 
tion of associations is 4/9 across both lines of sight), wc 
can estimate the likelihood of this result given the var- 
ious parameter sets. Assuming that each peak provides 
an independent test, and taking the probability of suc- 
cessfully finding a nearby quasar from Figure [T] (the long- 
dashed curves match the survey parameters reasonably 
well), then the probabilities to find exactly four associ- 
ations in the sample are (0.17%, 2.0%, 14%, 23%, 5.1%) 
for tQ = (0.3, 0.5, 1, 3, 10) x 10 7 yr, respectively, if we 
take ro = 45 comoving Mpc. Thus, with this single line 
of sight we can rule out at high confidence lifetimes of 
tQ ^ 5 x 10 6 yr (or even significant variability on those 
timescales) or tQ > 10 8 yr. 

However, taking a larger mean free path (ro = 200 
comoving Mpc) changes the constraints. In that case, 
the probabilities to find exactly four associations in the 
sample are (0.022%, 4.9%, 26%, 8.2%, 1.2%) for t Q = 
(0.1, 0.3, 1, 3, 10) x 10 7 yr, respectively We can still 
therefore rule out long lifetimes tQ > 10 8 yr by the lack 
of perfect associations, but we cannot place as strong 
constraints at the short end. Note that the constraints 
do not improve much more at large tQ, because the dis- 
tribution flattens out when tq exceeds unity for bright 
quasars. 

We also find probabilities of (0.11%, 8.5%, 9.5%, 0.5%) 
to find exactly four associations if f2 = 
(60°, 90°, 135°, 180°), respectively. 10 Thus here we 
can place strong constraints, ruling out any beaming 
scenarios with f2 < 90° as well as isotropic emission (at 
least with the infinite lifetime approximation); these 
are independent of the uncertainty in ro in our simple 
model where the nearest source provides a clear peak. 
Note that finite lifetimes and beaming both decrease 
the observed ratios, so combined models constrain small 
Q even more tightly (although isotropic emission with 
finite lifetimes is permitted). 

In this analysis we have used the low signal-to-noise 
Far Ultraviolet Spectroscopic Explorer (FUSE) He II 
spectra available at the time of the quasar searches 

it lies very close to one of the other two quasars, it actually adds 
no additional information. 

10 The small chance for isotropic emission (fl = 180°) accounts 
for quasars below the minimum luminosity limit. 
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(jWorseck fc Wisotzkil [200l IWorseck et all I2007D in or- 
der to al low a uniform sam ple across both spectra. Very 
recently, IShull et all (|2010l ) observed HE2347-4342 with 
the newly-installed Cosmic Origins Spectrograph (COS) 
on the Hubble Space Telescope. The remarkable improve- 
ment in the He II Lya forest spectrum clearly shows the 
peaks and troughs in the hardness ratio (their Fig. 7), 
and in principle this makes our test much easier. In this 
case it is very easy to identify peaks with a > 5; in fact 
the spectrum is so good that small-scale radiative trans- 
fer effects likely become important, since multiple peaks 
are often visible very close together. It is therefore some- 
what difficult to pick out the independent peaks in the 
spectrum, so a more detailed simulation is likely neces- 
sary in order to use this data most efficiently. 

For a simple estimate, we count the number of re- 
gions of extent r p ~ 20 comoving Mpc with peaks in- 
side them and use that as a proxy for the true num- 
ber of independent peaks. We count seven such regions 
in the range z = 2.4-2.72 (although several are very 
near our threshold); the spectrum cuts off below this 
rcdshift, and at higher values He II reionization may 
interfere with our simple model (iDixon fe Furlanettol 
l2009t iFurlanetto fc Dixonll2010t IShulT et al.ll2010D . Only 
one of the IWorseck et al.l (|2007ft associations sits in 
this range, so along this segment we have 1/7 peak- 
quasar associations. Using a = 5, we find the proba- 
bility of this result to be (4.4%, 24%, 35%, 8.0%, 0.43%) 
for tQ = (0.1, 0.3, 1, 3, 10) x 10 7 yr, respectively, and 
(39%, 27%, 3.2%, 0) or if ft = (60°, 90°, 135°, 180°), 
respectively. This spectrum is more tolerant of short life- 
times and tight beaming but less tolerant of long lifetimes 
and isotropic emission. As emphasized above, however, 
the data here are sufficiently good that simulations are 
clearly necessary to understand it. The difficulty of clas- 
sifying peaks and sources in this new data points to a 
statistical approach as most robust. 

For most of these quasars, the local minima in 77 co- 
incide in redshift with the quasars themselves, which 
implies that isotropic emission and relatively long life- 
times are good approximations. However, one source 
(at z = 2.69 toward HE 2347-4342) has a peak some- 
what in front of it according to the FUSE spectrum. 
IWorseck et al.l (|2007h estimate that t Q > 25 Myr from 
this coincidence. According to our models, the position 
in front of the quasar is suggestive of a finite lifetime as an 
explanation, but it could be explained equally well by a 
beaming model. On the other hand, the COS spectrum 
shows a peak at z = 2.69 as well as the one at lower 
rcdshifts, which would easily be consistent with longer 
lifetimes or isotropic emission. 

In any case, confusion from more distant sources as 
well as redshift errors can easily mimic this level of dis- 
placement in individual sources, so statistical samples 
are necessary for strong constraints based on the relative 
positioning of peaks and quasars. 

6. DISCUSSION 

We have examined how the transverse proximity ef- 
fect, observed through a combination of the He II and 
H I Lya forests, can help to constrain quasar properties. 
Using simple toy models, we showed how finite source 
lifetimes and beamed emission affect the statistical asso- 
ciation between peaks in the He II ionizing background 



and quasars. Both scenarios can substantially decrease 
the probability that the source causing a given peak is 
visible to a distant observer. Finite lifetimes break such 
associations when ctQ is less than the "proximity ra- 
dius" within which a quasar's radiation is more impor- 
tant than the accumulated background. Quasars with 
beamed emission may or may not be visible to both the 
Lya forest skewer and the observer for purely geometric 
reasons. As such, these models are less sensitive to the 
impact parameter of the quasar. 

We have shown that, even with the existing data (rela- 
tively small surveys around two lines of sight), this tech- 
nique suggests that tQ < 10 8 yr (with possibly some con- 
straints tQ > 3 x 10 6 yr if the older sample is more appro- 
priate for our toy model). Similarly, it appears to rule out 
either very small (< 60°) opening angles and isotropic 
emission (at least if an infinite lifetime is assumed). Of 
course, our simple toy models, which focus only on the 
single nearest source, are not sufficient to claim rigorous 
constraints. Improved numerical or Monte Carlo models 
can better test the importance of the accumulated back- 
ground of distant quasars, errors in quasar locations, and 
errors in peak detection in the He II forest. The im pres- 
sive COS He II spectrum from IShull et all (|2010ft illus- 
trates some of these difficulties, as one can clearly trace 
the peaks and troughs of the ionizing background, and it 
is difficult to determine which narrow peaks to associate 
with each source. Nevertheless, the constraints suggested 
by our toy model show that this is a very promising tech- 
nique. 

The mere presence or absence of a nearby quasar is 
also only the first bit of information from such surveys; 
the distribution of peak locations around quasars con- 
tains much more. Causality dictates that sources with 
finite lifetimes create peaks in front of (or at worst coin- 
cident with) the source, while beamed quasars can have 
their peaks behind the sources. The existing surveys have 
three peaks nearly coincident with their sources and one 
(possibly) leading; this provides weak evidence for finite 
lifetimes and/or anisotropy, but much better modeling is 
needed to interpret the data fully. 

In designing a survey, we have found that the bright- 
est neighbors are the most useful, because they have the 
largest proximity zones and produce the most obvious 
peaks. Faint quasars must be extremely close to the line 
of sight in order to create a peak; if they are so close, 
then the light travel time is small so they do not ef- 
ficiently constrain finite lifetimes. (The same is true of 
bright quasars very near to the line of sight; the strongest 
constraints on lifetimes will be provided by quasars at 
moderate distances from the Lya forest skewers.) 

We have argued that this method has two advan- 
tages over more traditional searches involving just the 
H I forest (fWorseck fc Wisotzkil 1200(1 ) . First, the trans- 
verse proximity zone is much larger in He II, because 
only (rare) quasars contribute to the high-energy ioniz- 
ing background and because the IGM attenuation length 
is several times smaller. Second, and more important, 
comparing the He II and H I Lya forests provides a di- 
rect measurement of the hardness of the radiation field, 
robust to variations in the underlying IGM density and 
temperature. This avoids a substantial bias in the H I 
proximity effect (jFaucher-Giguere et al. 1 120081) . 



Our tentative results from these toy models are consis- 
tent with the detection of the transverse proxim ity effect 
through metal lines by iGoncalves et all (|2008lh who es- 
timated Iq ~ 3 x 10 7 yr for two quasars. However , they 
may be inconsistent with lKirkman fc Tvtlerl ()2008f ) . who 
used a large sample of quasars with ~ Mpc separations 
to search the H I forest for the transverse proximity ef- 
fect. They found no evidence for enhanced transmis- 
sion in front of the quasars but decreased transmission 
behind them. They hypothesize that the increased gas 
density around the quasar hosts may cancel the expected 
increase in transmission in front of the quasars; in that 
case, the decreased transmission behind the foreground 
objects implies that the quasar light has not yet reached 
these regions, which in turn implies a short lifetime (or 
at least variability timescale) ~ 10 6 yr. Our constraints 
on variability on these very short timescales depend on 
the treatment of the data and so require more careful 
investigation. 

There are, in addition, many simplifications in our 
model. In addition to the inevitable measurement un- 
certainties in real experiments, we neglect the detailed 
structure of the forest and radiative transfer, which can 
induce small-scale features in the background. In the fi- 
nite lifetime case, the attenuation length of He II-ionizing 
photons is an important parameter, but one that is very 
difficult to measure directly. We have also assumed that 
optical surveys suffice to identify UV-bright quasars, but 
in reali ty there is substantial scatter in their far-UV prop- 
erties (|Telfer et al.1 l2002t IScott et all 12004 ) . which will 
degrade the correlation. 

Although we have focused on what the proximity ef- 
fect can reveal about quasar lifetimes and beaming, there 
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is more interesting physics to be gleaned. For example, 
seeing any proximity effect at all implies that a quasar 
can not be "flickering" more rapidly than the equilibra- 
tion time of the gas (~ 3 x I0 6 yr at T ~ 10 -14 s _1 , 
near the expected cosmic mean), or the mismatch from 
time delays wipes out the proximity effect. Because this 
timescale is much longer than that for H I, the resulting 
"flickering" timescale is actually significant compared to 
our expectations for quasar lifetimes. Thus the He II for- 
est can be used to measure long-term variability and per- 
haps even to constrain popular models like an exponen- 
tially decreasi ng quasar luminosity o r the more complex 
light curves of iHopkins et all ([2005). In either case, we 
might expect to find low-luminosity quasars (in their late 
stages of existence) correlated with surprisingly strong 
peaks in the hard ionizing background. This may be 
very interestin g in light of the rapid va riability timescales 
suggested by iKirkman fc Tvtlerl (|2008l) . 

Another possibility, in the context of finite quasar life- 
times, is to exploit the bare peaks by looking for a corre- 
lation between these and "post-quasar" galaxies, perhaps 
with a recently concluded burst of star formation or ev- 
idence for recent strong mechanical feedback from the 
quasar. Alternatively, in a beaming model, one would 
instead search for correlations with galaxies showing sig- 
natures of obscured quasars, such as rapid ongoing star 
formation. 
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